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Abstract. We present numerical results for the spin and thermal conductivity 
of one-dimensional (ID) quantum spin systems. We contrast the properties of 
integrable models such as the spin-1/2 XX Z chain against nonintegrable ones 
such as frustrated and dimerized chains. The thermal conductivity of the XXZ 
chain is ballistic at finite temperatures, while in the nonintegrable models, this 
quantity is argued to vanish. For the case of frustrated and dimerized chains, we 
discuss the frequency dependence of the transport coefficients. Finally, we give 
an overview over related theoretical work on intrinsic and extrinsic scattering 
mechanisms of quasi-lD spin systems. 

1 Introduction 

Quantum magnetism in ID is a successful example for a fruitful interplay between theory and 
experiment. On the one hand, many bulk materials exist that almost perfectly realize ID spin 
models (see |1I2) for a review) and on the other hand, powerful theoretical methods such as 
bosonization [3], the Bethe ansatz [4], or the density-matrix renormalization group method 
[5] are available. Often, excellent agreement between theoretical predictions and experiments 
has been found as far as ground-state properties, excitation spectra, thermodynamic or opti- 
cal properties are concerned (see |l|6j ). The understanding of transport properties is of great 
importance for the interpretation of transport or NMR measurements (see, e.g., Rcfs. [7|). Sub- 
stantial progress has been made in the past years, but the understanding is still incomplete, 
especially for systems involving many coupled degrees of freedom such as spins, orbitals, and 
phonons. Theoretically, the topic of transport in ID quantum magnets is challenging. First, 
several experiments demand for a more complete theoretical picture as will be outlined below, 
and second, transport theory often requires the computation of non-trivial correlation functions 
[5]. Third, transport is closely related to relaxation and non-equilibrium phenomena and thus 
connects to the rapidly evolving field of non-equilibrium physics of strongly correlated electron 
systems. 

In particular, the discovery of the colossal magnetic heat transport in spin ladder materials 
such as (Sr,Ca,La)i4Cu2404i, where the magnetic contribution to the total thermal conductivity 
k exceeds the phonon part substantially 9 10 11 , has sparked interest in transport properties of 
quasi-lD spin models. Often, a magnetic mean-free path is defined within a Boltzmann type of 
description [9|10j used to analyze the experimental data, which in the case of spin ladders can be 
of the order of several hundred lattice constants [10] . This observation - originally suggested to 
reflect ballistic transport properties of pure spin ladders - is not yet completely understood. 
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Several other spin-1/2 chain materials (see, e.g., [13|14j ) and 2D cuprate antiferromagncts 
[15|16j possess similar thermal transport properties, although typically the values measured for 
the magnetic contribution are much smaller. A very interesting aspect is the strong magnetic 
field dependence observed in some 2D |15|17j and ID materials |18|19j . Not all materials that 
exhibit a strong dependence of k on magnetic field are actually believed to have a significant 
contribution to k from magnetic excitations. Nevertheless, the possibility of tuning the thermal 
current through magnetic fields is appealing and may even allow to design functional devices 
such as spin valves [TS]. We refer the reader to a recent review and the article by C. Hess 
in this volume for more details on the experimental developments. 

Much theoretical work has focused on intrinsic transport properties of spin systems, address- 
ing intriguing questions such as the different transport properties of integrable as compared to 
nonintegrable ones. While in the remainder of this article, we restrict the discussion to the 
application of linear response theory - i.e., Kubo formulae - to spin and thermal transport of 
quasi-lD spin-1/2 systems as derived in Refs. [8121] . we note that alternative approaches such 
as master-equation techniques incorporating a modeling of heat baths have been pursued for 
quantum systems 22 23 24J. Moreover, while widely used, the derivation of Kubo formulae for 
heat transport may be questioned, as strictly speaking, no analogue to the voltage or magne- 
tization gradients driving electrical and spin currents exists in the case of thermal transport. 
We refer the reader to recent work on this issue [?4 25 26 27]. For brevity, we also concen- 
trate on spin-1/2 systems and refer to the literature for more details on Haldane systems [28] . 
Note, though, that due to similar low-energy properties [29], the transport behavior of gapped 
quantum systems such as spin ladders and spin-1 chains can be expected to be generic at low 
temperatures. Analogous questions, i.e., the properties of integrable vs nonintegrable systems, 
the validity of Fourier's law, and the modeling of heat baths are timely subjects in the study 
of transport of classical systems (see Ref. [30] for a review) . 

In linear response theory, ballistic transport is defined by the existence of a finite Drude 
weight D [31|32j . which is the zero- frequency contribution to the real part of the conductivity: 

Rck[o-](w) = D th[s] S(uj) + K[cr] leg ((j) , (1) 

where n denotes the thermal and a the spin conductivity. S(lo) is a <5-function and K[a] leg (uj) 
is assumed to be regular at u> = 0. Generally, the transport coefficients are computed from 
current-current correlation functions: 

or rco rf3 

k[<x]{w) = -^J dte^ +M + »j dT{j th[s] j th[s] (t + ir)). (2) 

Here and in all succeeding equations, r = for spin transport (labeled by 's') and r = 1 for 
thermal transport (labeled by 'th'). (3 = 1/T is the inverse temperature and (.) denotes the 
thermodynamic expectation value. A finite Drude weight implies a divergent dc conductivity. If 
D vanishes, then either a finite dc conductivity Ode = lim^^o c reg (a-0 can result, or, if <7 reg (oi;) 
exhibits an anomalous frequency dependence for uj — > 0, crd c may still diverge [33] . Note that 
here, we mainly consider finite temperatures, while the Drude weight was original introduced 
by Kohn to characterize a metal at T = |31|32j . 

Trivially, if the respective current operator commutes with the Hamiltonian, the Drude 
weight is finite at any temperature. It has long been known that the energy current operator 
of the spin-1/2 XX Z chain is a conserved quantity [M], but only later, a deeper connection 
between the existence of finite Drude weights at finite temperatures and the integrability of a 
model system has been made [35] . 

As a main objective of this paper, we wish to summarize recent theoretical progress, concen- 
trating on one-dimensional systems and their intrinsic spin and heat transport properties (see 
also Refs. [33136] for recent reviews). We will contrast the properties of integrable systems such 
as the spin-1/2 XX Z chain discussed in Sec. [2] against nonintegrable ones. As an example for 
the latter class of systems, we present numerical results for the spin and thermal conductivity 
of the frustrated and dimcrized spin-1/2 chain in Sec. [3] With respect to the experimental find- 
ings, obviously, both intrinsic as well as extrinsic scattering processes are of relevance. Recent 
theoretical results on extrinsic scattering channels are summarized in Sec. [4] 
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2 Transport properties of the XX Z chain 

We now turn to the nearest-neighbor spin-1/2 XX Z chain. The Hamiltonian is: 

1 



Hxxz — hi — J 



2 (S?Sr +1 + h.c.) + AS?S? +1 



(3) 



We set J = 1 in the following and periodic boundary conditions are imposed throughout this 
work. The current operators corresponding to the local energy density di = hi defined in Eq. ([3]) 
and local spin density di = Sf are obtained from the equations of continuity: 



N 



jth[s],m -jth[s],i = -*[-H>dth[s],i] => Jthfai = iy^\h-i,di] . (4) 



1=1 



It turns out that the energy current of the spin-1/2 XX Z chain is a nontrivial conserved 
quantity of this integrable model [35] Fl Hence in the case of thermal transport, Rck(cj) = 
D t hS(uj) for any exchange anisotropy A of this model. Although the spin current is not in 
general conserved in the case of the spin-1/2 XX Z chain, it has nevertheless been conjectured 
that D s should be finite at T > in the case of integrable models, but vanish in the case of 
nonintegrable ones [37I38| . While we will argue in Sec.[3]that this picture seems to be correct for 
the massive phases of nonintegrable ID spin models, some counterexamples have been proposed 
in the literature |39I40| . 



2.1 The thermal and the spin Drude weight at zero magnetic field 

As for the Drude weights of the spin-1/2 XX Z chain, the following picture has emerged: the 
thermal transport is ballistic for any exchange anisotropy and at all non-zero temperatures. Its 
dependence on T and A has been studied by means of Bethe-ansatz (BA) techniques |41l42j . 
exact diagonalization (ED) |12I43|44] . and with mean- field theory |43|45j . An example is shown 
in Fig. [lja), where we display ED data for the thermal Drude weight of the XX Z chain at 
A = 1 vs temperature [43] , in comparison with B A results from Ref. [44] . Note that numerically, 
the Drude weight can be computed from |35| 

D Ms] (T) = ^- £ e- E ^ T \(m\ Jth[s] \n)f. (5) 

Em=En 

Here, \n) and E n are eigenstates and -energies of H, respectively, and Z = ^2 n e ~ En l T denotes 
the partition function. Using system sizes as large as N = 20, the ED agrees with BA down to 
temperatures of T/ J ~ 0.25, which can be improved by employing extrapolation methods 1431 . 

Spin transport in the XX Z chain is a more involved problem as [H, j s ] ^ 0, and has been the 
objective of many studies |35l44l45l38l4fil47l48l49l5()l51l52l53l54l55l56l57l58l59lf)()l(ilj . While 
the spin Drude weight at T = is known exactly [62], for the massless regime \A\ < 1, where 
agreement exists that D S (T > 0) > 0, several BA calculations arrive at contradicting results 
for the temperature dependence |49|57|63] . The same holds for the question whether D s (T > 0) 
is finite or not at the SU(2) symmetric point A = 1. We refer the reader to Refs. [36157] for a 
discussion of this issue and the conceptual problems that BA approaches face. Recent numerical 
studies are consistent with D S (T > 0) > [44151161) . We illustrate this in Fig. QJb) and (c) 
where we show ED results for D S (T > 0) > as a function of temperature in Fig. [Tf b) and 
a finite-size scaling analysis of the leading residue C s = \iuiT^oo[T ■ D s ] in Fig. QJc). Within 
numerical precision and under the assumption that C s oc 1/N does not change at very large 
N, the extrapolation results in finite values for \A\ < 1. We refer to Ref. [44] for a detailed 
discussion of the finite-size scaling and to Refs. |44|50|56l for recent work on the massive, 
antifcrromagnctic regime A > 1. Numerical results for the ferromagnetic phase A < — 1 can be 
found in Ref. [43] . The frequency dependence of the regular part of the spin conductivity has 
numerically been studied in Refs. [38|48|56(64] . 



1 Trivial conserved quantities are, e.g., the total energy. 
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Fig. 1. Spin-1/2 XX Z chain, (a): Thermal Drude weight, ED (dashed lines: even N, dot-dashed 
lines: odd N; see [43]) vs BA results (solid line, [41]) at A = 1. (b): Spin Drude weight as a function 
of temperature and for several system sizes at A = 1. (c) Finite-size scaling at high temperatures: 
C s = limT^so [T ■ D s ] for A = 0, 0.5, 0.6, 1, 1.5 (from [44], with additional data for N = 20 at A = 1). 



2.2 Finite magnetic fields 

In the presence of magnetic fields Hz = &i > ^ ne ener Sy current jth and the spin current 

j s couple since particle-hole symmetry is broken and hence (jthjs) 7^ !35jEI This gives rise to 
magnetothermal effects, similar to the Seebeck effect of conduction electrons [8]. Transport in 
the XX Z chain in the presence of magnetic fields has been studied theoretically in [65|66l67|68j . 
focusing on the magnetothermal effect and the thermomagnetic power. To describe transport, 
now a 2 x 2 matrix description is necessary [8165166] : 



JA _ ( L n L 12 \ f Vh 
Ji ) \ L21 L 22 i I -VT 



(G) 



where Vh and VT are the gradient in field and temperature, respectively^ The conservation of 
jth is sufficient to show that all four Drude weights Dij corresponding to the transport coeffi- 
cients Lij are finite at all temperatures [35| . Then, assuming the condition of zero magnetization 
current flow, the thermal Drude weight -Kth is obtained as 

K th (h,T) = D 22 (h,T) ^ {k 'P T , , (7) 

where now the magnetothermal correction D\ 2 (h, T)/[T Dn(h, T)] contributes as well. Some of 
the main results of our work Rcf. [66] are (i) a reduction of K t h(h, T) due to the magnetothermal 
correction and (ii) expressions for the leading contributions to K t h(h,T) at low temperatures. 
The latter has been obtained from mean-field theory and bosonization. In the gapless phase of 
the XX Z chain (see, e.g., [5S] for details on the phase diagram), K t h{h,T) oc T, but at the 
quantum critical line separating the gapless from the ferromagnetic regime, we find K^ih, T) = 
AT 3 / 2 with A independent of the exchange anisotropy A. Note that Kt^h^T) has not been 
calculated yet by means of BA since D s (h,T) escapes an analytical treatment |57I67| . Mean- 
field theory, as outlined in Refs. [43|44|66j . proves useful as it provides a quantitatively good 
approximation to the thermal Drude weight both in zero [13] and finite magnetic fields [6T>] . 

In most cuprate based spin chain and ladder materials, exchange couplings are of the order 
of 1000K |l|2j hence little effects of a magnetic field on the thermal conductivity have been 
observed. In a recent experiment, the thermal conductivity of copper pyrazine dinitrate has been 
studied [19]. Since J ~ 10.3K [70], a significant field dependence is found at low temperatures. 
The analysis of the experimental data employing the mean-field theory description of Ref. [66] 
yields a constant mean-free path. The origin of this result especially at the quantum critical line 
remains to be elucidated. Magnetothermal effects do not seem to be present in this material. 



2 We invoke the notion of particle hole symmetry as ID spin models can be mapped onto spinless 
fermions with local interactions using the Jordan- Wigner transformation [8]. 

3 Note that Ji = j B and J2 = jth — hj B and therefore, D B — Dn. 
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Fig. 2. (a): Regular part of the thermal conductivity K reg (t<;) for a frustrated chain with a = 0.2 
[N — 18 sites; T/J =1,2 (circles, triangles)], (b): Finite-size scaling of the thermal Drude weight in 
the high-temperature limit: Cth = hm-r^oo [T 2 D t h] for a = 0.1, 0.2, 0.35, 1 (circles, squares, diamonds, 
triangles) for N = 8, 10, . . . , 20 sites (see also Ref. [44]). (c): Integrated spectral weight I t h for a = 0.2 
and T/J — 0.5,1,2. (d): Kreg(w) as a function of frequency u) for a = 1 [N = 18,20]. (e): Integrated 
spectral weight vs frequency at T/J — 0.3, 1,5 ((d), (e): reproduced from Ref. [72]). 

3 Transport properties of nonintegrable systems 

While originally conjectured to exhibit diffusive transport properties [37138171] , upon the exper- 
imental observation of large thermal conductivities in spin ladder, nonintegrable models have 
been discussed controversially in the literature |12|40|43|44|46|51|54;60|61|72|73|74|75rf6] . As of 
now, many studies point at a vanishing of both Drude weights [40 43 44 61172177178] in massive 
phases of nonintegrable models, including spin ladders. The masslcss regime of the frustrated 
chain remains a controversial issue [40|44j . Here we illustrate some numerical results for the 
finite-size scaling of the Drude weights [44] taking the example of the frustrated and dimerized 
chain and in particular, we also discuss the frequency dependence of the transport coefficients. 
The Hamiltonian of the dimerized and frustrated spin-1/2 chain is: 

N N 

2=1 1=1 

where a parameterizes the frustration and dimcrization is introduced through A; = 1(A) for an 
even(odd) site index I (A < 1). The current operators derive from Eq. ([4]) [43] , The regular part 
of k[ct](w) appearing in Eq. ([1]) can be written as: 

^HrcgH^I^ 1 '^ £ e~ E ^ T \(m\j th[s] \n)f6(u;-(E m ~E n )). (9) 
3.1 Frustrated chain 

As a result of preceding studies of the finite-size scaling of the thermal Drude weight [44], we 
concluded that no indications for a finite Drude weight are evident from the system sizes acces- 
sible by ED. This result is illustrated in Fig. HJb), where we show the leading coefficient of an 
expansion of the thermal Drude weight D t h(T) in powers of 1/T, i.e., C t h = limT^oo [7 l2 Z? t h (T)], 
as a function of system size (including new data for N = 20 as compared to Ref. [44]). The 
decrease of Cth with N is evident for all a as soon as the system size N becomes large enough 
(see Ref. [44] for details). The same picture arises for spin transport [44] . 

Let us mention the main features of Rc as found for the case of a = 1, i.e., in the 
massive regime, shown in Fig. [^d) [72] : (i) K reg (w) is a broad, featureless function extending 
up to frequencies ui/J < 4; (ii) at T/J = 1 and N = 20, the thermal Drude weight only gives 
a small contribution to the total weight of less than 3%. 

We now proceed by a discussion of the frequency dependence of the thermal conductivity of 
frustrated chains in the massless regime, i.e., a < 0.241 [79]. Our numerical results for n leg (u)) 
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Fig. 3. Dimerized chain with A = 0.1. (a): Regular part of the thermal conductivity as a function 
of frequency uj (N = 18 sites; T/J = 0.3,0.4,0.5; solid, dotted, dashed line), (b): Integrated weight 
Ith vs uj for T/J — 0.2,0.5, 1,2. Inset of (a): Integrated weight I t h(u) for uj/J — 0.5 as a function of 
temperature, (c): Regular part of the spin conductivity vs uj (N = 18 sites; T/J — 0.4,0.5, 1; solid, 
dotted, dashed line), (d): Integrated weight I s vs uj for T/J = 0.4, 0.5, 1, 2. Inset of (c): Enlarged view 
of the low-frequency region of panel (c). Vertical, dotted lines mark the position of the spin gap. 



and a = 0.2 are shown in Fig. [Ha) for N = 18 sites and T/J = 1,2. K Icg (uj) consists of a 
narrow peak centered around lo = 0, extending up to lo/ J < 0.05. This is reflected in the 
integrated spectral weight ith(w) = Jq dco' Rck(o/), depicted in Fig.^c). As this quantity also 
includes the contribution from the Drude weight, the figure reveals that the thermal Drudc 
weight on the system sizes considered here amounts to more than 50% of the total spectral 
weight. This observation is in stark contrast to the behavior of k(uj) in the massive regime on 
chains of a comparable length as summarized above. While on the one hand, the analysis of 
the finite-frequency properties of k for a > 0.241 supports the conclusion of a vanishing Drude 
weight, the question arises on the other hand whether the conclusion of D t h — > needs to be 
reconsidered in the massless regime. 



From the numerical data, it is difficult to establish a definitive conclusion about the ther- 
modynamic limit from small system sizes as far as current-current correlation functions are 
concerned. Still, no substantial differences are found between the massless and the massive 
regime concerning the finite-size scaling of the thermal Drude weight [44] . the common fea- 
ture being a monotonic decrease of Dth with system size N for N large enough. Moreover, 
bosonization studies on general grounds predict a vanishing Drude weight for thermal trans- 
port, irrespective of the presence of a gap [43|44|71(77] , In a recent work [78], Jung et al. have 
shown that, to first order in a, the commutator [H, j t h] between the Hamiltonian and the energy 
current operator vanishes, preserving the exact conservation of the energy current operator of 
the nearest-neighbor XXZ chain. While this feature explains the peculiar behavior of h(lu) 
for small system sizes as observed here, a vanishing of the Drude weight can still be expected 
in the thermodynamic limit where all terms of [H, j t h] in powers of a become relevant and 
cause a finite dc conductivity. The results of a quantum Monte-Carlo (QMC) study, however, 
seem to indicate that in massless phases of nonintegrable models, finite Drude weights may 
exist [40] . Note though that the interpretation of Monte-Carlo data at finite frequencies is quite 
involved as an analytic continuation from Matsubara to real frequencies needs to be performed 
|H9l4()l51j . 

In summary, the peculiar feature of k(uj) of frustrated chains in the massless regime, i.e. 
a < 0.241, as found for the system sizes accessible numerically, is that most spectral weight is 
found in the Drude weight, while the regular part consists of a narrow peak around lo = only. 



Will be inserted by the editor 



7 



3.2 The dimerized chain 

We next address finite-frequency transport properties of the dimerized chain (a = 0). In the 
following, we choose A = 0.1, i.e., we focus on the limit of strong dimerization A C 1. To 
first order in A, the dispersion relation of the elementary triplet excitation is described by 
£k/J = 1 + (A/2)cos(fc) [80], where k denotes the momentum. The spin gap G is quite large 
and roughly given by G / J = 0.95, while triplet-triplet interactions are suppressed by decreas- 
ing A. One may therefore on the one hand expect both the spin and heat conductivity to be 
small due to the large spin gap, but on the other hand, the transport properties should be well 
approximated by considering a weakly interacting gas of hardcore bosons [81) , which may, as a 
future project, allow for a comparison between numerical and analytical results. 

Our numerical results for the conductivities k(lo) and <j(lj) are presented in Figs.[3]Ja) and 
03c), respectively. The computations were performed for N = 18 sites, A = 0.1, and several 
finite temperatures as listed in the figure's caption. The distinctive features of both conductivi- 
ties visible in Figs.[3£a) and[3Jc) are: (i) Significant spectral weight is only found around u> = 
and in a high-frequency peak located around lo/J> 0.95, which corresponds to the spin gap. 
(ii) While the low frequency peak (including the Drudc weight ) contains a large fraction of the 
total weight in the case of the thermal conductivity, the spectral weight of the spin conductivity 
is mainly concentrated in the high-frequency peak. The latter is illustrated in Figs. Eljb) and 
Eljd), showing the integrated spectral weight I t h[s](w)//o, where Iq is the full spectral weight of 
K[er](u>). The low-frequency peak is present in a(ui) as well. The inset and the main panel of 
Fig. Etc) show that the low- frequency peak extends up to u>/ J ~ 0.1, which corresponds to the 
width of the one-triplet band. 

Furthermore, by integrating the low- frequency peak in /t reg (ii;) over w up to w/J « 0.5 
yielding I rcg (Lu/J = 0.5), we find that this quantity is independent of system size within nu- 
merical precision. I t h(uj/J = 0.5) is plotted in the inset of Fig. ]3fa) for N = 16 and N = 18 
sites. Hence, a significant redistribution of spectral weight as the system size increases is not 
expected. One further observes a maximum in J t h(w/J = 0.5) at roughly T/J ~ 0.35 and a 
1 /Independence at high temperatures. 

In summary, both models exhibit an intriguing behavior of the frequency dependence of 
both the spin and thermal conductivity that deserves further investigations. 

4 Extrinsic scattering 

As mentioned in the introduction, and as is evident form the phenomcnological analysis of 
experimental data for spin ladder |9|10|ll|83|82j as well as spin chain materials |13ll4ll84|85j . 
it is important to include external scattering processes to arrive at a realistic theory of thermal 
transport in quasi ID magnetic materials. For instance, doping with nonmagnetic impurities in 
(Sr,Ca,La)i4Cu2404i [82] - substitution of Zn for Cu - has been found to result in a suppression 
of the thermal conductivity linear in the Zn content. Mobile charge carriers effectively suppress 
the magnon thermal transport in spin ladder systems [83]. As heat transport via magnetic 
systems in a material requires the heat to be transferred from the lattice to the spin system, 
inevitably, spin-phonon scattering needs to be modeled by theory. 

First studies have addressed the thermal conductivity of spin-phonon coupled spin chains 
[77186187185] as well as spin ladders [55] • Some of these works [77189] start from effective field 
theories and describe transport within the Memory-matrix formalism [90191] by first identi- 
fying the slowest decaying modes, following the spirit of Rcf. [71]. These then determine the 
long-time behavior of current-current correlation functions. For the case of spin chains, an ex- 
ponentially large thermal conductivity K to tai oc exp(a(9/2T) is predicted [77], where is the 
Debye temperature. A peculiar result of the Boltzmann theory of Refs. [86187] is the constant 
spin thermal conductivity at high temperatures. For spin ladders, Rcf. 89J highlights the rel- 
evance of spin-phonon drag terms contributing to the total thermal conductivity, with a rich 
interplay of energy scales influencing the low-temperature behavior. A direct comparison of 
these results with experiments, however, needs to be done in future, in particular, as disorder 
may be of relevance in the structurally disordered spin ladder compounds (Sr,Ca,La)i4Cu2404i 
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[9|10|89j . Finally, note that spin phonon coupling has also been studied in the context of spin 
transport in the spin- 1/2 chain by means of QMC [92] . 

Note that via the Jordan- Wigner transformation Heisenberg type of models can be mapped 
onto spinless fermions [8], the transport properties of which have extensively been studied in 
the context of localization [S3]. We just mention an incomplete list of recent, closely related 
works addressing Heisenberg chains [94 95 96 97 98 99], spin ladders [73], or effective low-energy 
models ]86|87|100] . Interestingly, some works seem to indicate that the dc spin conductivity may 
be finite for interacting systems in the case of off-diagonal disorder [95196] , Also, even if the dc 
spin conductivity vanishes, the same is not necessarily true for thermal transport as energy can 
still be transfered over a weak link [100] . 

Finally, only results from a mean-field theory are available for the thermal conductivity of 
doped spin ladders in the literature [101] , Transport properties of ID t-J and Hubbard models 
have widely been investigated (see, e.g., Ref. [33] for an overview), and it is beyond the scope 
of this work to discuss the charge and spin transport of these systems. Their thermal transport 
properties have, however, not been studied sufficiently [102] . Note that the Hubbard model, 
being integrable, is expected to exhibit ballistic thermal transport, which also holds for the 
supersymmctric point of the t-J model [35] , 

5 Summary 

We may conclude that the intrinsic thermal transport properties of the spin-1/2 XX Z chain 
in zero and finite longitudinal fields are well understood. The spin transport of this model still 
poses some challenges to theorists, such as an analytical calculation of the spin Drude weight 
of the spin-1/2 Heisenberg chain. As for nonintcgrablc systems and within linear response the- 
ory, it seems that generically, ballistic transport in the sense of finite Drude weights is not 
realized. Rather, the relevant information is encoded in the frequency dependence of the con- 
ductivities. The challenge to computational scientists is to devise algorithms that can simulate 
low temperature regimes. Analytical approaches face the problem that effective field-theories 
of nonintegrable models are typically integrable, with diverging transport coefficients. Hence, 
the definition of a low-energy theory that describes transport accurately is a nontrivial task. 

Promising results with respect to the interpretation of experiments have been obtained from 
first studies incorporating phonons or disorder, but a consistent picture has not emerged yet. 

Highly interesting and potentially new physics is expected from both experiments and novel 
theoretical methods such as the time-dependent density matrix renormalization group method 
[103] that investigate transport and relaxation of strongly-correlated electron systems away 
from equilibrium. 
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